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B 

C 
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i(t) 


K(|,i|r a ,t b ,t c ) 


s 

T 
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Definition 

beam A, upstream beam 

beam B, downstream beam 

a constant (C = 41^/otD) 

diameter of laser beam; photodetector 

absolute time averaged current supplied to photo- 
detector with knife-edge removed 

ac-electr ical signal output of photodetector 

similarity function 

cross-correlation function or cross-correlogram 
of beam A 

sensitivity of laser schlieren beam 

integration time for computation of cross-correlation 
time 


<U> the most probable speed of disturbances averaged over 

the beam separation distance | 

x,y,z coordinates in direction of flow, perpendicular to 

flow direction (in horizontal plane), and perpendicular 
to flow direction (in vertical direction), respectively 

QI Glads tone-Dale constant 

A beam deflection vector 


| beam separation distance 

p air density 

T time delay of downstream signal 

0 angle of rotation of A during transit time of disturbance 

from beam A to B 


iv 



DEFINITIONS OF SYMBOLS (Cont'd) 


Symbol 


<w > 

“P 

( )' 
T ) 
LI 
C) 

< > 


Definition 

angle measured from x-axis to beam deflection vector, A 

the most probable preferred angular velocity of the 
disturbances 

fluctuating component 

time average 

vector 

unit vector 

expectation value, or most probable value 


Subscripts 
x,y, z 
a,b 
c 

1,2 

m 


component in direction of x,y,z-axis, respectively 
beam A,B 

refers to angle between beam B and the Z-axis 
case A, no rotation; case B, rotation 
most probable; maximum 
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A DIRECT MEASUREMENT OF THE MOST PROBABLE 
PREFERRED ANGULAR VELOCITY OF TURBULENT STRUCTURES 
BY OPTICAL CORRELATION OF LASER SCHLIEREN SIGNALS 

SUMMARY 


A method is introduced which provides the means for separating 
the translational and rotational motion of turbulent' s tructures . Simple 
two-dimensional models are used to relate the skewness of cross-correlo- 
grams computed from laser schlieren signals to the rotation of flow 
disturbances. The method, referred to herein as the method of "forced 
similarity," is discussed with respect to application to the turbulent 
free shear layer of an axisymmetric supersonic jet. Experimental results 
show that the shape of the cross-correlogram in the neighborhood of the 
"peak" is strongly influenced by rotational motion, and therefore, it 
becomes necessary to account for the effect in order to determine the 
correct statistical properties of the turbulence. 


I . INTRODUCTION 


This document contains the author's preliminary thoughts on the 
possibility of making direct measurements of the most probable preferred 
angular velocity at "localized" regions within turbulent flows by optical 
correlation of laser schlieren signals [2]. Although some experimental 
results are presented which satisfy certain necessary conditions required 
of the proposed theory, sufficient data verifying the theory do not pre- 
sently exist. However, plans have been made to obtain these data in 
Marshall Space Flight Center's Cold Flow Thermal and Acoustic Jet 
Facility. 

During the feasibility test presented in reference 2, there were 
unexplained variations in the shape of the cross -correlograms computed 
from laser schlieren signals retrieved from the supersonic turbulent 
boundary layer on a thin plate. In April 1968, near the conclusion of 
this test, the connection between the skewness of the cross -correlograms 
and the angular rotation of the disturbances was considered. 



Preliminary experimental results eventually led to a method for 
separating rotational and translational contributions to the computed 
cross -correlograms . This method is referred to herein as the method 
of "forced similarity." The "forced similarity" condition provides the 
means for direct measurement of the most probable preferred 1 angular 
velocity of the turbulent structures 2 at localized regions inside turbu- 
lent flows. The potential extension of the method to yield the probabil- 
ity distribution of the angular velocity and measurements of vorticity 
are very interesting, and may result as a natural consequence of the 
development of the laser schlieren system. 


II. THE LASER SCHLIEREN OPTICAL REMOTE SENSING SYSTEM 

For simplicity, let us consider the model of a laser schlieren 
optical remote sensing system shown in figure l 3 where two laser beams 
of light are directed through the test section of a wind tunnel. The 
parallel beams are perpendicular to the flow and separated by a dis- 
tance |. The plane formed by the beams is such that flow disturbances 
passing through the upstream beam (beam A) at time t, also pass through 
the downstream beam (beam B) at a later time t + t. A knife-edge is 
positioned perpendicular to each of the beams such that 50 percent of 
the light is prevented from reaching the respective photodetectors when 
the beams are undisturbed. 

Reference 2 shows that a disturbance characterized by a local gra- 
dient of the index of refraction will deflect each of these beams as it 
passes through them. Each deflection is proportional to the component 
of the gradient which is perpendicular to the path of the beam. When a 
deflection occurs, the amount of light reaching the photodetector is 
changed. The photodetector converts this change into an ac-electr ical 


x The word "preferred" is used since the statistical process does not 
yield the angular velocity components based upon statistical averages 
of absolute values. 

throughout this report reference is made to "turbulent structures" and 
"disturbances" in order to distinguish between the motion of the density 
gradient and the motion of fluid particles which represent an "eddy" at 
a particular instant. The laser schlieren signals are produced by the 
fluctuating density gradient component perpendicular to the path of 
the beam [3] . 

3 This system is described in detail in reference 2. 
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signal. Therefore, the output signal of the photodetector is related 
to the beam deflection which is caused by a local gradient character- 
izing a particular flow disturbance passing through the laser beam. 

Because the refractive index of air is proportional to the density 
(to a good approximation), the output signal of the photodetector is 
proportional to the change in the component of the density gradient 
which is perpendicular to the beam and to the knife-edge. 


i(t) 


- cJs^eLM' 

ox 


avg. 


( 1 ) 


In equation (1), s is the sensitivity of the beam, a is the Gladstone- 
Dale constant [3], and i(t) is the output signal of the photodetector. 
Also, the knife-edge was assumed to be perpendicular to the x-direction 
as shown in figure 2. 

Figure 2, a view taken along the centerline of the laser beam 
(either beam A or B) from laser toward detector, shows the knife-edge, 
the eye of the photodiode, and the laser beam cross section. The beam 
is .shown in the undisturbed position and in a deflected position 
described by the beam deflection vector, A. 

In the following let A(t) and 0(t) represent the magnitude and 
direction of A(t), respectively. Further, assume that the eye of the 
photodiode has a large area of constant sensitivity compared to the 
cross sectional area of the laser beam, and that no laser light falls 
outside of this area of constant sensitivity. If the intensity across 
the laser beam is assumed to be constant over the beam cross-section 
and that the knife-edge is straight, the relationship between the time 
histories of the beam deflection vector and the photodetector ac-elec- 
trical output is simplified considerably. 

Because the area of sensitivity of the eye of the photodiode 
monitoring the fluctuating laser light is constant, the component of 
the beam deflection vector which is parallel to the knife-edge does not 
contribute to the output signal of the photodetector. The component of 
A which is perpendicular to the knife-edge will determine the amount of 
change in light reaching the photodiode. If the beam deflection is 
small compared to the diameter of the laser beam, the output signal, 
i(t), can be expressed conveniently as a function of the component A 
which is perpendicular to the knife-edge: 


i(t) 



A (t) 

X 


( 2 ) 
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UNDISTURBED BEAM 


Figure 2 


A - BEAM DEFLECTION VECTOR 


A Schematic of A Cross Section Taken Across 
One Beam of A Laser Schlieren System. 





where the first group of terms in brackets on the right-hand side of 
equation (2) is the change in the output signal per unit change in A x 
and is assumed to be a constant, C. 


i(t) = C • A x (t) 


(3) 


For this case, the photodetector output is a linear function of the com- 
ponent of A which is normal to the knife-edge. This result is not 
unrealistic and can be approximated very well in practice. 

Since 


A (t) = A(t) • cos 0(t), (4) 

equation (3) can be expressed as 

i (i) = C • A(t) • cos 0(t) . (5) 


Thus, we see from equation (5) that i(t) is a function of the direction 
as well as the magnitude of the beam deflection vector. 

For the laser schlieren system shown in figure 1, the signals from 
the photodetectors monitoring the respective beams are 


i a (t) = C a ' A a (t) * cos 9 a (t) ( 6 ) 

and 

i b ( fc ) = c b ’ ^(t) - cos e b (t). (7) 

Let us assume that the turbulence is two-dimensional; i.e., the 
time-averaged statistical properties do not vary along the beams. 
Although this assumption will restrict the analysis, it provides a 
better model for the purpose of describing the fundamental concept, 
the prime purpose here. It will soon become evident that relaxing 
this assumption will not change the fundamental relationship between 
the shape of the correlograms and the angular velocity of the density 
gradient vector component which is normal to the respective beams. 
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by 


The cross -correlation of the signals i a (t) with i^Ct + x) is given 


T 

R(|,t) = lim ; f i (t) • iu ( t+ T> dt. (8) 

T J a b 

o 

Let us substitute equations (6) and (7) into (8) and assume that the 
integration time, T, is large enough so that 


R(I,t) 


C C, 
a b 


A a (t) 


(t+x) • cos 0 a (t) * cos 0^(t+x) dt. 


(9) 


Figure 3 shows the knife-edge orientation for the cross-correlation of 
equation (9). This figure also shows the beam deflection vector of beam 
A at time t when a particular disturbance is passing through the beam, 
and, the A of beam B when the same disturbance is passing through beam 
B at a later time t + t. When this disturbance passes through beam A, 
the. beam is deflected in a direction 0 a , and by a magnitude of A a - 
This magnitude of the deflection vector is determined by the magnitude 
of the density gradient component normal to the beam characterizing the 
disturbance. The direction of the deflection is due to the particular 
orientation of the same density gradient component. If we make the 
realistic assumption that the disturbance rotates as it travels down- 
stream, its orientation will not necessarily be the same when it passes 
through beam B as it was when it passed through beam A. This is repre- 
sented in figure 3 where the beam deflection vector caused by a particu- 
lar disturbance has rotated during its transit from beam A to beam B. 

It follows that 


A a (t)x Aj^ (t+x) 


A a (t) . Ah (t+x) • sin 0(t+x). 


( 10 ) 


This equation defines the angle 0 through which the disturbance rotated 
during transit from beam A to B. Thus, 


0(t+-r) = e a (t) - e b (t+x). 


(ii) 
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Solving equation (11) for 9 b (t+r) and substituting into equation (9) 
gives 


C a C b f 

R(£,t) = / A a (t) • A fe (t+T> • cos 0 a (t) . cos[0 a (t) - 0(t+T)]dt. 

o 

( 12 ) 

To examine the effect "which this angular rotation has upon the 
cross-correlation function R(|,t), let us consider a simplified two- 
dimensional flow model where it is assumed that the disturbances are 
random and produce statistically stationary signals. The cross-corre- 
lation function (equation (12)) will be studied for two cases; (a) where 
the disturbances do not rotate and (b) where each disturbance rotates 
through the same angle, 0 , during transit from beam A to beam B. Case 
(a) will be discussed first. 

Figure 4 shows four knife-edge arrangements where the view is taken 
along the laser beams similar to that of figure 3. The differences 
between these four figures is the orientation of the downstream knife- 
edge. If we assume that the disturbances do not rotate, then 


0(t+T) = 0, 

and from equation (11), we see that 


(13) 


e b (t+x) = e a (t). 


(14) 


The signal from beam A will be 


i a (*-) = C a • A a (t) • cos e a (t) 


(15) 


and is the same for all four knife-edge arrangements. 

The signal from beam B, for figures 4a,b,c, and d are 

i bl( t+T ) = C b ' A t,( t+T ) ' cos 0 b ( t+T ) (16) 
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1 b2 ( ' t+T ' ) 

- c b ' ' 

s in 

e b (t+T) 

(17) 

l b 3 (t+T) 

= -C b • A b (t+ T ) 

• cos 

e b (t+x) 

(18) 

i b4 (t+T) 

- ‘ C b ' A b< t+ ^ 

• s in 

6 b ( t+T ) > 

(19) 


respectively. 

When we substitute equation (14) into equation (16), (17), (18) 
and (19), the cross -correlation functions for the respective knife- 
edge arrangements become 


Ri(S,t) = 


c c, 

a b 


Aft) • A, (t+x) • cos 0 (t) • cos 0 (t) dt (20) 
a b a a 


R 2 (!,t) = 


c c, 

a b 


o 

o 


A (t) ' A b (t+ T ) • cos 9 a (t) • sin 0 g (t) dt (21) 


R3(6,t) = 


c c, 

a b 


A a (t) • A b (t+ T ) 


cos 0 (t) • cos 0 (t) dt (22) 
a a 


R-id.'O = - 


c c - 

a b 


A (t) • A, (t+x) cos 0 (t) • sin 0 (t) dt. (23) 
a d a a 


The cross-correlograms for these arrangements are shown in figures 5a, 
5b, 5c, and 5d, respectively. An example of a cross-correlogram, repre- 
sented by equation (20), would be similar to that shown in figure 6a. 

A cross-correlogram similar to that represented by equation (22) is 
shown in figure 6b. 
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CROSS-CORRELATION 


TIME DELAY 


( b) 



TIME DELAY 

Figure 6(a). Typical Cross -Correlogram for Knife-Edge 
Arrangement of Figure 4(a) (Equation 20) 
(b). Typical Cross-Correlogram for Knife-Edge 
Arrangement of Figure 4(c) (Equation 22) 



For case B, it is assumed that each disturbance rotates through 
the same angle 0 during transit from beam A to B; then, 


0(t+x) = 0- (24) 

Expanding equation (12) and making the substitution from' equation (24), 
we get 


R(i,x) = 


c c, 

a b 


A (t) • A, (t+r) • cos 2 0 (t) • cos 0 dt 
a b a 


+ 


-C C, 
a b 


A (t) • A, (t+x) • cos 0 (t) • sin Q (t) • sin 0 dt. 

3 b 3 3 

(25) 


Since 0 is not a function of time, cos 0 and sin 0 can be taken outside 
of the integrals on the right-hand side of equation (25), respectively. 


R(|,t) = [cos 0] 


C c, 

a b 


A (t) • A, (t+x) • cos 2 0 (t) dt 
3D 3 


+ [sin 0] 


C C,q 
a b 


o 

T 


A & (t) • A b (t+x) • cos 9 a (t) • sin 9 a (t) 


dt. 

(26) 


In equation (26), the integrals represent the cross -correlations derived 
for case A. Substitution of equations (20) and (21) into (26) yields 


R(|,x) = [cos 0] Rid.x) + [sin 0] R 2 (e,x). 


(27) 
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Equation (27) relates the cross -correlogram of case B to those of case 
A. We see from equation (27) that 


R(S,t) = Ri(|,t), 

for 0=0 

(28) 

R(i,x) = r 2 (6,t), 

for 0 = jt/2 

(29) 

R(£,t) = -Ri(I.t), 

for 0 = it 

(30) 

R(I,t) = -R 2 (|,t), 

for 0(3/2) rc 

(31) 


as should be expected. 

In equation (27) the unknown quantity is the angle 0. A direct 
measurement of 0 for the idealized model in case B can be made by com- 
bining case A with case B. This procedure is referred to here as the 
method of "forced similarity" and is as follows; 

(1) Compute the cross-correlogram R(£,t) for zero beam separation 
with the knife-edges as shown in figure 2 ((■ = 0). Since 0=0 for 
| = 0 , 


R(0, t) = R 1 (0, t) . 


(32) 


(2) Separate the two beams by moving beam B downstream a known 
distance | from beam A. The position and knife-edge angle, ij/ a , of 
beam A are not changed (\|r a = 0) . 

(3) Rotate the downstream knife-edge, B, to the angle which 
produces the same "shape" of the cross-correlogram as was computed for 
zero beam separation. For this simplified flow model, the angle _ 
corresponding to the maximum "degree of similarity" between the cross- 
correlograms is equal to the angle 0 through which the disturbances 
rotate during transit from beam A to beam B. The "forced similarity" 
condition is 


R(!,t - VV = • R(0,t) (33) 

where K(|, x, ij/ b ) is "optimum" when ty b = i|r bm for a particular 

(4) The most probable transit time x m of the disturbances is 
determined by the time delay on the cross-correlogram for separated 
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beams corresponding to the similar position at zero time delay on the 
cross-correlogram for zero beam separation. 

(5) The most probable speed of the disturbances is 


<U > = | t . 

s m 


(34) 


(6) The most probable angular velocity is 


<w > = 0T = \k x . 
— p 7 m T bm m 


(35) 


These models have been used to relate the experimentally observed skew- 
ness of cross-correlograms , computed from laser schlieren signals, to 
the rotation of the turbulent structures. 


III. PRACTICAL APPLICATION OF THE METHOD OF FORCED SIMILARITY 


In the previous section a method of "forced similarity" was intro- 
duced which provides the means for measuring the most probable preferred 
angular velocity of the simplified two-dimensional disturbances. As a 
consequence of the proposed relation between the skewness of the cross- 
correlogram and the rotational motion of the disturbances, the "forced 
similarity" condition, equation (33), must be imposed in order to deter- 
mine the most probable transit time, T m , as well as all other statis- 
tical properties which are computed from the shape of the cross-correlogram 
in the neighborhood of x m . Imposing this condition analytically involves 
normalization of the cross-correlogram and optimization of the similarity 
function, K(£,T,\|q>). These details will not be discussed here. Rather, 
we will proceed to the practical application. 

Figure 7 shows a schematic of the free jet shear layer of a super- 
sonic axisymmetric air jet. Two laser beams are directed through the 
flow perpendicular to one another and in such a manner that the plane 
formed by the beams is perpendicular to the centerline of the jet. The 
horizontal beam (beam A) passes through the center of the jet. The 
vertical beam (beam B) passes through the shear layer and intersects the 
horizontal beam as shown in figure 7. 

In reference 4 it is theoretically shown that the statistical cross- 
correlation of the signals from two such beams should result in a cross- 
correlogram representative of the components of the signals which are 
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Figure 7. Cross-Beam Arrangement for Axisymmetric Jet 


caused by disturbances passing through the localized region about the 
beam intersection. That is, flow disturbances passing through the beams 
that are not common to both beams do not contribute significantly to 
the cross-correlogram (the integration time being sufficiently long). 

In reference 2 the experimental verification of this "cross-beam" concept 
was successfully achieved in a supersonic (M = 2.0) turbulent wake of a 
thin flat plate. 

Referring again to the axisymmetric jet, the "forced similarity" 
condition may be satisfied for retrieving all three components of the 
most probable preferred angular velocity <qu>, of the disturbances 
passing through a localized r.pgion of the shear layer. We may assume 
that the mean statistical properties of the turbulence are also axisym- 
metric, and thereby eliminate one component of <Wp>. However, this is 
not necessary and will not be done. 

The experimental procedure is as follows; 

(1) The beams are directed through the turbulence as shown in 
figure 7 with zero beam intersection (| = 0). The knife-edges are 
positioned as shown (ijr a = 0 and ^ = 0) , 

(2) The cross-correlogram R(0 ,t) is computed. 

(3) The vertical beam (beam B) is moved downstream a distance | 

as shown in figure 8, and the cross-correlogram R(|,t) is computed which 
satisfies the "forced similarity" condition 


r (I,t - V^am’W = K ^> T >t am >W * R(0 ,t) (36) 

for axisymmetric flows. K(£ , T, t am , jr bm ) is determined to be the "weakness" 
function necessary to satisfy equation (36) in the neighborhood of t . 
Experimentally, this involves rotation of both of the knife-edges. m ' 

When this step is completed the condition should be satisfied for this 
particular flow. 

(4) If necessary, the next step would be to rotate the downstream 
beam in the plane perpendicular to the centerline of the jet. This is 
done in increments and step (3) is repeated for each increment until all 
three angles \Jr a , tJtj-, , and \)r c have been determined which satisfy the three- 
dimensional similarity condition 

R (i>T " Vtam’^bm^cm) = K(lj T>t am^bm’W * R(0 ’ t) < 37 > 


where t)/ c is the angle of rotation of the vertical beam. 
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(5) The most probable transit time of the disturbances is deter- 
mined by evaluating the similarity condition (37) for t = 0; 


R(|,0 - T ,i|r , \|/, , ik ) = 

m Y am Y bm T cm 


R(0,0) 




’em' 


(38) 


Thus, T m is equal to the time delay on the cross -correlogram R(|,T,i|/ am , 
^bm’^cm) corresponding to the similar position at zero time delay on the 
cross -correlogram for zero beam separation, 

(6) The components of the most probable preferred angular velocity 

are 


<w x> 

T cm 

T 

(39) 


m 


<(*) > = 

“\lr 

Y am 

(40) 

y 

T m 

<10 > = 

!bm 

(41) 

z 

T 

m 


and the most probable vector is 

<to > - i <co > + j <w > + k <io >. (42) 

— r\ -v *■' v 


(7) The most probable speed of transit is 
< U > = . 


(43) 


Figures 9 and 10 show the first successful cross-beam measurements 
made in the supersonic free shear layer of an axisymmetric jet. The 
exit Mach number of the nozzle was 2.5, and the expansion was optimum. 
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CROSS - CORRELATION CROSS-CORRELATION 




Figure 9. Cross -Correlogram Computed from Laser Schlieren Signals 
Retrieved from Free Shear Layer of Axisymmetric Jet 
(M e =2.5). 

(a) Positive Time Delay Range 

(b) Negative Time Delay Range 
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CROSS-CORRELATION (volts 



Cross -Correlogram Computed for Same Case as That of Figure 9 
Except Beams are Separated by 1.0 Inch as Shown in Figure 8 
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In figures 9a and 9b, the positive and negative time delay ranges 
of the cross -correlogram are shown. The beam geometry is shown in 
figure 7 (| = 0) . It can be seen that the correlogram is not symmetric 
about the origin (t = 0), but that there is a dominant peak at (t = 0). 
This cross -correlogram corresponds to R(0 ,t) in equation (36). 

Figure 10 shows the cross -correlogram for a beam separation of one 
inch (| = 1") with the knife-edges in the same positions as in figure 9 
(i|/ a = 0, = 0) . The fact that there is a considerable difference in 

the shapes of these correlograms satisfies one necessary condition 
required by the theory of "forced similarity." 

Figure 11 shows four cross-correlations computed from signals 
retrieved with a laser schlieren system using parallel beams as pre- 
viously described. The data" were retrieved from the supersonic 
(M = 2.0) turbulent boundary layer on the thin plate shown in figure 12. 
The beams were separated by approximately 1.5 inches in the direction of 
flow and were approximately 1/8 inch above the surface of the plate. 
These cross-correlograms were computed with the downstream knife-edge 
orientations shown and were the first experimental attempt to relate 
the shape of the correlogram, in the neighborhood of the most probable 
transit time, t , to the orientation of the knife-edges relative to the 
flow and to each other. 

These data show that the relative orientation of the knife-edges to 
each other does indeed influence the shape of the cross-correlogram. 
Therefore, it is reasonable to consider that the cross-correlogram is 
influenced by a change in the relative orientation of the disturbances 
to the knife-edges as they translate from one to the other. 


IV. CONCLUSIONS 


Based upon the theoretical and experimental results which have been 
presented, the following conclusions are; 

(1) The relative orientation of the knife-edges to each other con- 
trols the shape of the cross-correlogram in the neighborhood of the time 
delay, x m , corresponding to the most probable transit time of disturb- 
ances between the laser beams. 


These measurements were made in MSFC' s Bisonic Wind Tunnel on June 3, 
1969 [2], 
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Figure 12. Shadowgraph of Flow Field Generated by the Thin Plate Model 
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(2) The method of "forced similarity" introduced here provides 
the means for separating the translational and rotational contributions 
to the shape of the cross-correlogram. However, more experimental 
results are needed for verification, 

(3) The "forced similarity" condition is based upon the assump- 
tion that the shape of the cross-correlogram computed with beams 
separated is "similar" to that computed with zero beam separation when 
the contribution due to rotation has been eliminated (the effect of 
decay accepted) . 

(4) Because of the rotational influences, significant error can 
result in the flow properties which are computed from the shape of the 
cross-correlogram unless the influence due to rotation is considered. 

(5) The condition of "forced similarity" may be applied by 
rotation of either or both beams about the axes of a reference 
coordinate system in combination with rotation of the knife-edges. 

The particular combination used is primarily dependent upon the flow. 
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